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ABSTRACT 

We present the preliminary results of Monte Carlo simulations aimed to investigate 
the effects of realistic dust extinction (absorption + scattering) on the colours of high-z 
galaxies. In this paper, we concentrate on the case of spheroidal galaxies, and we obtain 
attenuation curves in the range 0.1-1.2/im for different dust spatial distributions, and 
for a range of values of the dust optical depth, geometrical thickness and inclination 
of the dust disk. We find that the resultant curves are strongly dependent on the dust 
geometrical distribution and optical depth. A serendipitous finding is that the strength 
of the 2200 A absorption feature depends not only on the optical depth, but also on 
the dust geometrical distribution. As a first application, we test our results on two 
high-z galaxies with extremely red colours (HRIO and HR14) in order to infer clues 
on their ages, dust content and dust spatial distribution. We confirm that HRIO must 
be a very dusty galaxy, and we suggest that its stellar component should be strongly 
embedded in the dust in order to reproduce the observed extremely red colours. 

Key vifords: 

dust: extinction ~ radiative transfer - scattering - galaxies: general - galaxies: evolu- 
tion - galaxies: individual (HRIO, HR14) 



1 INTRODUCTION 

Internal dust might affect seriously the observed proper- 
ties of distant objects, and observations have indeed found 
evidence of dust in high-z objects. Several active galax- 
ies with z >2 have been detected at submillimeter wave- 
lengths, suggesting the presence of large amounts of dust 
(Mdust ~ 10*~^ M0; see Hughes 1996 for a recent review). 
The reddening of the background quasars and their metal 
abundances suggest the presence of dust in damped Lya 
absorption systems (Pettini et al. 1994; Pei & Fall 1996). 
Also, the UV polarization properties of high-2 radio galax- 
ies can be explained in term of dust scattering, suggesting a 
significant amount of dust in their ISM (Cimatti 1996 and 
references therein). A substantial amount of dust is also ex- 
pected in evolutionary models of spheroidal galaxies at high- 
z (Franceschini et al. 1994; Mazzei & De Zotti 1996 and 
references therein). 

Although dust is likely to be present in most high-2: 
systems, no information is available neither on its spatial 
distribution, nor about how realistic dust extinction can af- 
fect our view of distant galaxies. In fact, extinction is usu- 
ally treated in a simplistic way, neglecting the contribution 
of scattering, and assuming naive spatial distributions (e.g. 



uniform foreground screens or infinite slabs). Only recently, 
more realistic models have been developed (Kylafis & Bah- 
call 1987; Bruzual et al. 1988; Witt et al. 1992; Byun et 
al. 1994; Wise & Silva 1996; Bianchi, Ferrara & Giovanardi 
1996, hereafter BFG). Although these models have shown 
that dust scattering plays an important role by reducing 
the effects of dust absorption, no systematic studies of the 
effects of the extinction on the observed colours of high-z 
galaxies have been performed. Understanding these effects 
is crucial in cosmology and galaxy evolution studies. 

For instance, a relevant problem arises in the age es- 
timates of high-z galaxies. When deep spectroscopy is not 
available, no information on the stellar continuum and ab- 
sorption features are obtainable, and the age estimates are 
based solely on the fitting of the broad-band photometric 
Spectral Energy Distributions (SEDs) with synthetic stellar 
population spectra. In a simplistic scenario where the ex- 
tinction is entirely due to absorption (the screen model), it 
is well known that the colours of a reddened young galaxy 
can mimic those of an unreddened old galaxy, producing a 
degeneracy of the age and dust extinction. However, what 
happens to the integrated colours of a galaxy when dust 
extinction is treated more realistically is poorly known. 

Witt et al. (1992) suggested that the effects of realis- 
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tic dust extinction on the SEDs of high- 2; galaxies are not 
strong. However, their claim was based on a single object 
(the radio galaxy B2 0902+34 at 2; ~ 3.4), and generalized 
to the whole population of high-z galaxies. In addition, the 
continuum SED of B2 0902+34 turned out to bo completely 
different from that used by Witt et al. (1992) because of a 
strong contamination of the /f-band flux by the redshifted 
[OIII]AA4959+5007 (Eisenhardt & Dickinson 1992). 

Instead, recent studies (Franceschini et al. 1994) have 
emphasized the role of dust in the early stages of evolution 
of galEixies. The present optical surveys for high- 2; galaxies 
are designed to select objects with strong emission lines, 
such as Lya, (sec Pritchct 1994 for a review; Thompson et 
al. 1995), or galaxies with a fiat continuum spectrum and 
with a sharp Lyman-break (Steidel et al. 1996). However, 
because of those selection criteria, these surveys would miss 
the putative population of high- 2; galaxies obscured by dust 
extinction. It is notable that the strongest evolutionary ef- 
fects due to dust extinction are not expected for the galac- 
tic disks, or in disk dominated systems (Mazzei et al. 1992), 
but rather in early-type galaxies which, under appropriate 
circumstances, might experience a prolonged opaque phase 
(Mazzei & Do Zotti 1996). The whole issue of the effects of 
dust on the SEDs of high- 2; galaxies still remains an impor- 
tant, open question, especially in view of the future ISO and 
sub-mm data. 

Extinction effects become particularly relevant in the 
case of high-z galaxies with very red colours, also called ex- 
tremely red objects (EROs) (McCarthy et al. 1992; Hu & 
Ridgway 1994). In fact, these galaxies may be distant and 
old ellipticals, and could provide crucial clues on the first 
epoch of galaxy formation, Ho and go, once their ages are 
derived accurately. However, they could also be very dusty 
high-2 galEixies whose intrinsic colours are strongly reddened 
by dust extinction. The main problem is that the degener- 
acy of the age and the dust extinction becomes maximized 
in this class of objects because of the similarity between 
the colours of a genuinely old galaxy and those induced by 
foreground screen reddening. The main questions are then : 
how much realistic dust extinction can mimic the colours of 
an old galaxy ? If those galaxies are dusty, is it possible to 
learn something about the spatial distribution of the dust 
by modeling their SEDs ? 

Motivated by the general lack of information on the ef- 
fects of dust extinction on the SEDs of high-z galaxies, we 
have started an extensive study aimed at investigating the 
relevance of these effects. In this paper, we focus on the case 
of spheroidal galaxies, present the first results, discuss an ap- 
plication to the case of the extremely red galaxies, and show 
a serendipitous result about the strength of the 2200 A dust 
absorption feature in external galaxies. More general results 
and applications will be presented in a forthcoming paper. 



2 MONTE CARLO SIMULATIONS 

As a first attempt to study the effects of dust extinction, 
we investigate the case of spheroidal galaxies. The stellar 
density profile pstar{r) is modelled as a ,]affc bulge (sec 
BEG), which reproduces the surface brightness r^^'* pro- 
file characteristic of elliptical and bulge systems (de Vau- 
couleurs 1959). We adopt an effective radius re =4 kpc, rep- 



resentative of nearby elliptical and bulge systems (Binney 
& Tremaine 1987), and the distribution is truncated at a 
radius rmax = 5re = 20 kpc (see BEG for more details). 

Dust extinction has been treated by using Monte Carlo 
simulations of realistic radiative transfer in dusty galaxies 
(i.e. considering the extinction as the combination of ab- 
sorption and scattering). The details of the model are not 
repeated here and can be found in BEG. However, it is im- 
portant to recall here that we assume Galactic dust: in our 
model the grains are supposed to be spherical and to have a 
size distribution given by the MRN model (Mathis, Rumpl 
& Nordsiek 1977), n{a) oc a~^'^, where a is the grain ra- 
dius. We consider three materials: astronomical silicates, || 
graphite, and ± graphite. The numerical silicates/graphite 
ratio is 1:1.12, with 1/3 of the graphite having optical prop- 
erties measured parallel (||), and 2/3 perpendicular (+) to 
the c-axis. The lower and upper limits of the distribution 
are a_ = 0.005 /im and a+ = 0.25 fim, irrespectively of the 
material. Very small grains and PAH have not been included 
due to the large uncertainties both in their size distribution 
and optical constants. The dielectric constants adopted are 
the ones given by Draine & Lee (1984), recently extended in 
the far UV and X-rays by Martin & Rouleau (1991). All the 
relevant optical properties (absorption and scattering cross 
section, albedo) have been calculated using Mie formulae. 

We tested three different spatial distributions for the 
dust within the galaxies : (A) a homogeneous disk with con- 
stant dust density pdust and with radius = rmax = 20 kpc 
equal to that of the stellar distribution. We considered three 
different values of the geometrical half-thickness of the disk, 
Zd = 70, 250, 500 pc, and three different inclinations to the 
line of sight, i = 10°, 45°, 90° (edge-on); (B) a homogeneous 
spherical distribution with radius — rmax = 20 kpc equal 
to that of the stellar bulge, where the dust is intermixed with 
the stars and has constant density {pduat{r) =constant); (C) 
a spherical distribution where the dust is intermixed with 
the stars, but its density follows the same radial distribu- 
tion of the stellar component (i.e. pdust{r) = Pstar{r)). 

Model A can be considered representative of a 
spheroidal galaxy where the dust is located in a disk, sim- 
ilarly to what is actually observed in a rmmbcr of nearby 
ellipticals (Goudfroij 1996, and references therein). On the 
other hand, in models B and C, the dust is interspersed 
with the stars. In particular, model B, although simplistic, 
can mimic a very dusty spheroidal galaxy where the stars are 
strongly embedded in the dust, and it is somehow similar to 
the classic foreground screen model, i.e. where the dust has a 
large geometrical covering factor, and where dust absorption 
dominates. Finally, model C can be considered representa- 
tive of a spheroidal galaxy where the dust is distributed 
more realistically, following the same gravitational potential 
as the stellar component. We are aware that high-« galaxies 
may have different distributions of dust, and that models A, 
B and C have simplistic geometrical configurations. How- 
ever, these models allow us to test on the observed SEDs of 
high-z galaxies a broad variety of circumstances, from very 
dusty cases (model B) to cases where the dust is located 
only in a disk (model A), and have the capability to treat 
dust extinction realistically, i.e. taking into account both 
absorption and scattering. 

In addition to the above geometrical parameters, each 
model is characterized by the value of the optical depth in 
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the V-band in the rest frame of the galaxy, tv, along a line 
of sight passing through the centre of the galaxy for models 
A and B (perpendicular to the plane, in the disk case). In 
the case of model C, due to the central cusp of the r^^^ law, 
it is preferable to define rv for a line of sight at projected 
distance 1.16 re from the centre (see BFG). The simulations 
were carried out for a range of tv from 0.1 to 10, and for 28 
wavelengths in the range 0.1-1.2/im. 



3 RESULTS 

3.1 Attenuation curves 

From the Monte Carlo simulations we obtain the attenua- 
tion curves ^Ia (defined as the ratio of the observed, dust ex- 
tinctcd, total intensity, lobsW to the intrinsic, unextincted, 
one of the gala:xy, /o(A)) as a function of wavelength for 
different values of tv, Zd, and i (see Fig. 1). 

As expected, in the case of the disk (model A) with 
i = 90°, we find an almost flat A\ for any value of tv and 
Zd (i.e. no effects on the colours). For inclinations i < 90°, we 
find instead a strong dependence of on tv, whereas the 
depencence on Zd, and i is much weaker. A representative 
case, Zd=250 pc and i — 45°, is shown in Fig. 1: the shape of 
the curve changes significantly from tv -^0. 1-1.0 to tv >1.0, 
and the depth of the 2200 A feature varies considerably, 
reaching a maximum for tv ~1, becoming weaker for larger 
Tv, and almost disappearing for tv > 4. Regarding the case 
of the disk, it is important to note that we found that an 
exponential distribution of the dust within the disk does not 
produce any relevant difference with respect to the case of 
a disk with a constant dust density. 

For the bulge with constant dust density (model B), also 
shown in Fig. 1, the slope steepens considerably for tv > 1, 
and the depth of the 2200 A feature shows remarkable vari- 
ations as a function of tv- Moreover, at a given tv, model 
B produces attenuation curves much redder than in model 
C. 

In sunnnary, for a given tv there are appreciable differ- 
ences between the three models in terms of the slope of the 
.4a and strength of the 2200 A feature. A detailed analysis 
of these results will be presented in a forthcoming paper. 

3.2 The 2200 A feature and its absence in 
external galaxies 

The absence or weakness of the 2200 A absorption feature 

has been recently debated extensively (Calzetti et al. 1994 
and references therein; see also Calzetti et al. 1995). The 
problem arises when the 2200 A feature is expected, but not 
observed in galaxies which are known to be dusty from their 
emission line ratios and UV-optical continuum properties. 
The absence of the feature has been usually ascribed either 
to a different dust chemical composition or to scattering 
effects which compensate for the dust absorption. 

Our study shows that a peculiar chemical composi- 
tion of the dust grains, or the role of scattering may not 
be the ordy causes to explain the lack or weakness of the 
2200 A feature in external galaxies. Figure 1 clearly shows 
how the strength (measured by its equivalent width) of the 
2200 A feature is a strong function of both the optical depth 



Figure 1. A representative set of attenuation curves obtained for 
the three models : model A (solid), B (sliort-dashed), and C (long- 
dashed). The curves are sliown for tv=0. 1,0. 5,1. 0,2.0,4.0,10.0, 
with Tv increasing from the top to the bottom of the figure. 



and of the geometry. The feature shows strong variations as 
a function of the optical depth within each of the models 
A, B and C. For instance, in model A (disk case) its equiv- 
alent width varies remarkably depending on tv, reaching a 
maximum value of ^ 500 A at tv ~ 1 and decreasing to 
~ 65 A at TV ~ 6. In addition, strong variations are also 
evident, for a fixed tv, as a. function of the geometry of the 
dust extinction. For instance, for = 1, the 2200 A fea- 
ture displays large differences among the models A, B and 
C. In this regard, it is important to note that, at least for the 
disk distribution (model A) , the feature is strongly weakened 
or suppressed because the emission from opaque galaxies is 
mostly contributed by stars outside the dust disk and hence 
unaffected by absorption. 

The strong variations of the feature as a function of 
the geometric configuration suggest that the compensat- 
ing effect of scattering may not be the only way to reduce 
the strength of the feature, and that the geometry may 
have a dominant influence. It is important to recall here 
that the main contribution to the bump comes from small 
(a < 0.1/um) graphite grains; for the dust composition here 
adopted, the weighted average radius at this wavelength is 
in fact ~ 0.05 nm. Correspondingly, the albedo is rather low, 
(w) ~ 0.2, and the scattering phase function rather isotropic, 
{g) ~ 0.2. Thus, it seems unlikely that scattering is the main, 
or the only, origin of the weakness of the 2200 A feature 
compared to the pure absorption case. However, a more de- 
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tailed analysis is required to answer this question unambigu- 
ously. 

Although our results suggest that the strength of the 
2200 A feature may be a strong function of the geometrical 
effects, they do not rule out the possibility that external 
gala^Kies may have chemical grain properties different from 
those of our Galaxy. In fact, it is important to recall that 
the 2200 A feature is weak in the LMC, and absent in the 
SMC (Prevot et al. 1984). The issue of the 2200 A feature 
will be investigated in details in our forthcoming paper. 



4 APPLICATION TO EXTREMELY RED 
GALAXIES 

As a first application, we investigate the effects of dust 
extinction on the SEDs of two extremely red galaxies at 
high-z. Extremely red objects, EROs, (i.e. R — K > 6) arc 
found both in random sky fields or in AGNs fields (Elston 
et al. 1988; McCarthy, Persson & West 1992; Hu & Ridg- 
way 1994). Since their colours may be strongly affected by 
dust extinction, these galaxies represent an interesting ap- 
plication of our models, and suitable objects to investigate 
at what level the SEDs of dusty younger galaxies can mimic 
those of old ellipticals once we treat dust extinction in a re- 
alistic way, and to constrain the spatial distribution of the 
dust. 

We have selected from the literature two extremely red 
galaxies, HRIO and HR14, for which the best optical-IR pho- 
tometric SEDs are available (Hu & Ridgway 1994; Graham 
& Dey 1996). Our study is more concentrated on HRIO, the 
only ERO with spectroscopically determined redshift. 

We attempted to fit the SEDs of HRIO and HR14 with 
synthetic stellar population spectra, and applying the at- 
tenuation curves derived by our models in order to take 
into account the possible dust extinction occurring in these 
galaxies. In our experiments we made use of the synthetic, 
unextincted, stellar population spectra of Bruzual & Char- 
lot (1993), with instantaneous burst of star formation, solar 
metallicity, and Salpeter IMF. 

Figure 2 shows our results for HRIO, a galaxy at z=1.44 
with a colour index I — K =6.5. Since its quoted B-band 
flux is not statistically significant (~ 2.2ct level), we adopt 
a 3(7 upper limit. We also correct for the ^^20% H-band flux 
contamination from the Ha line. Graham & Dey (1996), 
adopting a foreground screen dust absorption model, found 
that the SED of HRIO cannot be reproduced without dust 
extinction. Our independent analysis confirms that dust ex- 
tinction must play an important role in HRIO, and even 
very old (i.e. very red) stellar populations (up to 7 Gyr) 
cannot reproduce its SED without extinction. However our 
modeling allows us to go a step further, and to derive clues 
on the spatial distribution of the dust in HRIO. In fact, we 
find that model A (dust disk) cannot fit the data because 
the observed colours are always too red compared to those 
produced by this model for any tv, age, Zd, and i. Model 
C allows a better fit to the SED, but the reddening is still 
not sufficient to match the I — K colour. We find that only 
model B gives an acceptable fit to the SED (age ~ 1 Gyr and 
Ty=6). The SED can be fit also by a 1.0 Gyr old popula- 
tion reddened by a dust foreground screen with tv=1.54 (or 
Eb-v=0-55, adopting the Galactic relation Ay ~ 3Eb-v)- 



Figure 2. The rcst-framc UV-optical SED of HRIO. Flux data 
points arc from Graham & Dey (1996). The curves are synthetic 
stellar population spectra cxtinctcd by dust. The observed SED 
cannot be fit without extinction, but only model B is capable to 
produce enough reddening to match the colours of the galaxy. 

If we take these results at face value, we infer that the dust 
in HRIO is not distributed in a disk, and may be interspersed 
with the stars in a way that dust absorption dominates, re- 
sembling the effects of a simple screen model, i.e. the stars 
may be strongly embedded in the dust. This also shows how 
such modeling may provide new clues on the internal dust 
distribution in unresolved, distant galaxies. The Ha emis- 
sion line properties observed by Graham & Dey (1996) sug- 
gests that HRIO may be a very dusty galax;y with star for- 
mation and/or AGN activity. Adopting the model B with 
TV — G and a 1 Gyr old stellar population, we estimate 
that the rest-frame dereddened K-band absolute magnitude 
is MK„,t=-27.5 {Ho = 50 km s"^ Mpc~^, 90=0.5), about 
9 times brighter than Af^ of the luminosity function of the 
local field galaxies (Mobasher, Sharpies & Ellis 1993). 

Figure 3 illustrates the results relative to HR14, a 
galaxy with / - K =6.2 (Hu & Ridgway 1994). The case 
of this galaxy is more ambiguous because its spectroscopic 
redshift is not known. Hu & Ridgway (1994) estimated a 
photometric redshift Zphot='2.3, with a ±3(t acceptable fit 
range 1.8< z <3.0. Although its redshift is uncertain, we 
have decided to investigate the case of this galaxy in three 
different cases : z=1.8, 2.3 and 3.0, which represent the ac- 
ceptable range of photometric rcdshifts. For 2=1.8, we find 
that its SED can be reproduced successfully without any 
dust extinction and with very old stellar populations (>4 
Gyr). However, if dust extinction is taken into account, the 
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Figure 3. The observed SED of HR14. Flux data points are from 
Hu & Ridgway (1994), and the case of 2=2.3 is shown. The curves 
are synthetic stellar population spectra which provide acceptable 
fits to the observed SED. The SED of HR14 can be also fit with 
an unextincted old stellaj: population spectrum. 



HRIO, the only ERO with a useful set of data to attempt to 
answer this question. 

First of all, it is necessary to check whether the pho- 
tometry of HRIO may be contaminated by strong emission 
linos, such as [OIIIJAA4959+5007 in the ,J-band. Let us sup- 
pose that HRIO is a heavily obscured quEisar. Adopting a 
typical ratio [OHI]/Ha <0.25 observed in Seyfertl /quasars 
(Netzer 1990), and using the observed Ha flux (Graliam & 
Dey 1996), we find that the J-band flux contamination is 
small (<8%). We attempted to fit the SED of HRIO with 
a reddened quasar spectrum and a foreground dust absorp- 
tion screen model. For this purpose, a point-like source (the 
quasar nucleus) embedded in a dusty medium, the screen 
model is a good approximation for the dust extinction. For 
the quasar spectrum, we used the average spectrum of radio- 
quiet quasars of Cristiani & Vio (1990) {Fn oc v~'^'^) be- 
cause HRIO is a radio-quiet object. We find that the SED of 
HRIO cannot be reproduced with a reddened quasar spec- 
trum for a wide range of Eb-v- We also tried to use flatter 
{Fi, oc v~^''^) and steeper {Fi, oc v^^'"^) quasar spectra in 
order to take into account the dispersion of the slope of the 
average quasar spectrum, but we always found that the SED 
cannot be reproduced by a reddened quasar spectrum. 

Therefore, although the Ha-|-[NII] line luminosity 
(3x 10'*^/ij~gg erg s~^; Graham & Dey 1996) may be com- 
parable to that of Seyfert 1 and quasars, our analysis of the 
SED suggests that the spectrum of HRIO is not dominated 
by AGN radiation. However, the large FWHM of the line 
(7000± 3000 km s"^) seems to favour the presence of an 
AGN, or to indicate that part of the radiation comes from 
an active nucleus with broad lines. Deeper observations are 
needed to establish the value of the FWHM, which is only 
marginally significant. 



SED can be fit within the errors by models A, B and C, and 
with a wide range of stellar population ages and tv (for in- 
stance, model A: tv ^ 3.0, age 1.6 Gyr; model B: tv ^ 2.0, 
age 1.6 Gyr; model C : - 2.0, age 2 Gyr). For z=2.i, the 
results are shown in Figure 3. The main result does not differ 
much from the case of z=l.d, : shifting HR14 to z=2.?, has 
the effect of requiring younger ages for all cases compared 
to 2=1.8. It is also interesting to note that the foreground 
screen (absorption only) requires Tv ~ 0.84. For z=Z.Q, we 
find results almost identical to those of 2=2.3, except for the 
Tv=0 case, where the SED is best fit with a slightly younger 
age (2.5 Gyr instead of 3.0 Gyr). 

As a general result of this study, we note how realistic 
dust extinction (i.e. absorption -|- scattering) requires larger 
dust optical depths compared to the case of the pure ab- 
sorption screen model (see also Witt et al. 1992). This is 
evident for example in Figure 1, where in order to produce 
the same amount of reddening, it is required Ss_v=0.55 
(i.e. TV ^ 1.5), and tv ~ 6 for the screen model and model 
B respectively. 

4.1 Can HRIO be a dusty quasar ? 

Hu & Ridgway (1994) found that the surface density of 
EROs is comparable to that of quasars. This result raises 
the question of whether EROs axe dusty quasars and AGN 
missed by optical surveys. We concentrate this discussion on 



5 SUMMARY AND CONCLUSIONS 

We have presented the preliminary results of a set of simu- 
lations aimed to study the effects of realistic dust extinction 
on the colours of high-z galaxies. As a first attempt, we have 
considered three different spatial distributions of the dust in 
spheroidal galaxies, adopting Galactic dust grain properties. 
A serendipitous finding is that the 2200 A feature is gener- 
ally much weaker than in the foreground absorption screen 
model. The weakening of the feature is particularly evident 
in the case of a disk of dust in a spheroidal galaxy, where the 
feature can be almost absent. We find that the equivalent 
width of the 2200 A feature depends on the dust distri- 
bution, geometry and optical depth. We suggest that this 
effect is mostly a geometrical one, and that it might be one 
explanation of the failure of the detection of the 2200 A fear 
ture in external dusty galaxies. However, other causes of the 
non-detection of 2200 A feature, such as different chemical 
composition of the dust grains, cannot be ruled out by the 
present results. 

As a first application of our extinction curves, we con- 
sider the case of two extremely red galaxies at high redshift, 
HRIO and HR14, in order to constrain the ages of their 
stellar populations, the amount of dust, and the dust geo- 
metrical distribution. Our results confirm that HRIO must 
be a very dusty galaxy, and allow us to constrain the spar- 
tial distribution of the dust in this object. However, we find 
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that a strong degeneracy between the age of the stellar pop- 
ulation and the dust extinction is present in both galaxies. 
In the case of HRIO, we find that its colours can be repro- 
duced by a 1 Gyr old stellar population and tv ~6 (model 
B). Its dereddenod rest-frarnc K-band absolute magnitude, 
Mifrest=-27.5, is just beyond the highest luminosities of the 
local field galaxies. 



For HR14, although the redshift is poorly constrained, 
we find that its SED can be reproduced without extinction 
and with an old stellar population (2.5-4.0 Gyr) over the 
wide range 1.8< z <3.0. This makes HR14 a good candidate 
of a mature galaxy at z >1.8. However, if dust extinction is 
present, tiie age can be lowered down to 0.8 Gyr for rv 
3 and for a broad range of dust spatial distributions. It is 
important to note that the main results on HR14 are not 
strongly dependent on z, if 1.8< z <3.0. If HR14 would 
turn out to be at the same redshift of HRIO, we would find 
results very similar to those obtained for HRIO because of 
their similar I — K' colours. We also explore the possibility 
that HRIO is a heavily obscured quasar, but we find that its 
SED cannot be reproduced by a reddened quasar spectrum. 
However, the present data cannot exclude that part of the 
radiation in HRIO is due to an AGN. 



Since realistic dust extinction requires larger optical 
depths compared to screen models (see also Witt et al. 1992), 
it is important to estimate the total amount of dust required 
by our models. For standard values of the average grain size 
(a ~ 0.1/im) and density (5 ~ 2 g cm~^) the total dust mass 
in these galaxies is in the range Mdust ~ SaR^axTv = 10®~® 
Mq for TV = 1 — 10, respectively. In particular, taking 
into full account the dust composition, our models have 
Mdust ~ krv 10* M0 , where k depends on the model ge- 
ometry { kA ~ 1.21, fcs = 0.81, kc ~ 0.35). These masses 
are comparable with those derived for several active high-z 
galaxies by submillimetric observations (Hughes 1996). 

The existence of a population of extremely red galax- 
ies, whose colours can be explained only in terms of strong 
dust extinction, may be relevant to our understanding of the 
galaxy evolution. In fact, a dusty phase is predicted to occur 
during the early evolution of spheroidal galaxies (Frances- 
chini et al. 1994; Mazzei & De Zotti 1996 and references 
therein; Zepf & Silk 1996). If this is correct, a substantial 
population of dusty galaxies is expected to be present in 
the early universe. Those galaxies would bo missed by the 
current optical selection criteria used to find high- 2; galaxies 
(see Pritchet 1994,PASP,106,1052; Steidel et al. 1996). How- 
ever they would be selected by deep surveys in the sub-mm 
spectral range, where the redshifted dust thermal emission 
should peak (Blain & Longair 1993). It is tempting to spec- 
ulate whether the population of very red galaxies recently 
found by deep optical and near-IR imaging (Hu & Ridgway 
1994; Knopp & Chambers 1997; Hogg et al. 1997) may be 
somehow related to this putative family of high- 2: galaxies 
hidden by dust obscuration. Future multiwavelength obser- 
vations of the extremely red galaxies will provide new clues 
on their nature and role in the galaxy evolution picture. 
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